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ABSTRACT 



The spacings between bound-state levels of the Schrodinger equation 
with the same principal quantum number N but orbital angular mo- 
menta i differing by unity are found to be nearly equal for a wide range 
of power potentials V = Xr u , with Em ~ F(u, N) — G(u, N)£. Semi- 
classical approximations are in accord with this behavior. The result 
is applied to estimates of masses for quarkonium levels which have not 
yet been observed, including the 2P cc states and the ID bb states. 

I. INTRODUCTION 

The properties of heavy quarkonium (cc and bb) levels have provided useful 
insights into the nature of the strong interactions. At short distances these interac- 
tions appear to be characterized by a Coulombic potential associated with asymp- 
totically free single-gluon exchange Jjj], while at large distances the interquark force 
is consistent with a constant, corresponding to a separation energy increasing lin- 
early with distance [0. An effective power potential V ~ r u , with v close to zero, 
provides a useful interpolation between these two regimes for cc and bb levels when 
the interquark separation ranges between about 0.1 and 1 fm 0. 

Charmonium (cc) levels have been identified up to the fourth or fifth S-wave, 
the first P-wave, and the first D-wave excitation. (We are ignoring fine structure 
and hyperfine structure for the moment.) Six S-wave levels and two P-wave levels 
have been found in the upsilon (bb) family. Many additional levels are expected but 
have not yet been seen. There are some reasons of current interest for predicting 
their positions in a relatively model-independent manner. For example: 

(1) It has been suggested || ||] that, depending on its exact mass, the second 
charmonium P-wave level (which we shall denote as ^ C (2P)) could play a role in 
the hadronic production of the ip'(2S) state. 

(2) The first D-wave level of the bb system, which we shall call T(LD), may be 
accessible to experiments of improved sensitivity at the CLEO detector, both in 
direct production of the 3 Di states via the e + e~ channel and through electromag- 
netic transitions from the T(3S) state 0. 

Explicit calculations in nonrelativistic quark models of the masses of the x' c (2P) 
and T(1D) states have tended to have very small spreads || [TJ. Any interquark 
potential which reproduces the known quarkonium levels is well-enough determined 
to leave little room for variation in predictions of these levels. However, in the 
course of examining these predictions within the context of power potentials V = 
\r", with Xu > 0, we were struck by a curious feature: The spacings between levels 
with the same principal quantum number but orbital angular momenta differing 
by unity are nearly identical for a wide range of values of v. 
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The principal quantum number N is that which labels nonrelativistic Coulomb 
energy levels through the Balmer formula E^g oc — 1/iV 2 . It is related to the radial 
quantum number n r and the orbital angular momentum £ via N = n r + i (n r = 
1,2,3, . . .), where n r — 1 = n corresponds to the number of nodes of the radial 
wave function between r = and r = oo. Our result amounts to a formula for 
levels E N g linear in £ for a given N: 

E m = F{v,N)-G{v,N)l (1) 

with G(v,N) > for v > — 1 in accord with a much more general result which 
specifies the order of levels for fixed N as a function of I [§]. This general result is 
that G is positive for any potential whose Laplacian is positive everywhere outside 
the origin. We investigate this general case first in Section II, turning to the special 
case of power potentials in Sec. III. We apply the results to specific quarkonium 
cases in Sec. IV, while Sec. V summarizes. Some proofs of identities are given in 
three appendices. 

II. POTENTIALS WITH POSITIVE LAPLACIAN 

We consider the Schrodinger equation in appropriate dimensionless units (h = 
2fi = 1, where \x is the reduced mass), for a spherically symmetric potential V(r): 

[-A + V(r)]y m {r) = E m y m {r) . (2) 

Potentials with positive Laplacian AV > play a special role in quarkonium 
physics. First, all potentials used to fit the bb and cc spectrum have positive Lapla- 
cian (except for spurious pathologies when one expands the potential in powers of 
the strong coupling constant at a given scale). It is in fact very natural to take 
potentials with positive Laplacian because this property is a kind of expression of 
asymptotic freedom: If we say that the force between a quark and an antiquark is 
F = —a(r)/r 2 , asymptotic freedom requires that a(r) be an increasing function 
of r. Writing F = —dV/dr, saying that a(r) increases is equivalent to 

dr dr 

which means that V has a positive Laplacian. If this is the case we have 0, for a 
given multiplet, 

En,£-1 > Enj > En,£+i (4) 

for N > 3, t > 1. 

Very soon after the discovery of this property, Baumgartner || proved the 
following theorem: 
Define 

Rn,i=~^ p • (5) 
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Then 

if V = — (1/r) + Xv, with Av > 0, A sufficiently small. Notice that the right-hand 
side of (|6|) is always less than unity. So (|6|) holds for perturbations of a Coulomb 
potential. It is not known if this inequality also holds outside the perturbative 
regime. Our belief is that it does hold, but the proof must be rather complex. The 
perturbative proof is based on the equation 

£ati(E N/ -i - E N/ ) -{£+ l)a 2 {E N/ - E N/+l ) 
2/4-1 r°° 

= —^1 AV(r)u 2 N/ dr (7) 

where Un,£ is a pure Coulomb wave function, and 

Aa 2 e = (£ + I)' 2 - N- 2 . (8) 

For AV > 0, the inequality (H follows from (|7|) and (f|). In fact, inequality @ 
would also hold for a potential with purely negative Laplacian, since both ([?]) 
and (^) are reversed. The proof of (0) will be given in Appendix A. What is 
important is that this equation shows that Baumgartner's result || cannot be 
improved, because it is saturated by a potential v whose Laplacian has its support 
concentrated at the zeroes of u n ^. Such a potential is easy to construct: It is 
made of piecewise shifted Coulomb potentials. The resulting overall potential is 
not concave, but can be made so by a correction of order A 2 (A being the coefficient 
of v), for example by joining pairs of Coulomb segments by linear pieces tangent 
to each. 

Even if one believes that (|6]) holds for non-perturbative situations, it is a rather 
disappointing result for practical applications. For instance, using the "particle 
physics" spectroscopic notation, we get 

E 2 p — E\d 16 

while, as we shall see in the next section, this ratio can be very close to unity. On 
the other hand no improvement is possible. Even the very simple-minded potential 
— (1/r) + Ar gives 

E-ss — E 2 p _^ 1 
E 2 p — Eip, 2 

for A — > 0, as shown in Appendix A, and we have checked numerically that for 
finite A this limit is approached smoothly. Therefore, to make progress, we should 
use a more restricted class of potential, the power potentials, which are known to 
give excellent fits of the heavy quark-antiquark spectra H. 
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III. THE CASE OF POWER POTENTIALS 



If V(r) = e[y)r v ', where e{y) is the sign function and — 1 < v < 2, one finds 
that Rn,i is always very close to unity, to summarize briefly this section. The evi- 
dence comes from indications from semi-classical formulae valid for large quantum 
numbers, analytical investigations in the neighborhood of v = — 1 and v = 2, and 
explicit numerical calculations of energy levels for small quantum numbers, shown 
in Fig. 1. For the purpose of this figure, in order to obtain a smooth limit as v — ► 0, 
we have represented the Schrodinger equation as 



-(l/2)A+(r"-l)/i/]*Ur) 



(11) 



and levels are labeled by the radial quantum number n r = N — £ and the orbital 
angular momentum £ = (0, 1, 2, 3, . . .) =(S, P, D, F, . . .). 

For power potentials, semi-classical formulae are expected to hold for iV — > oo. 
However there are several regimes depending on whether £ stays small compared 
to N or, on the contrary, N — £ stays small. In the latter case two of us QTOJ have 
shown that depends asymptotically only on the combination N — £(1 — ^l +2 ) ■ 
In the former case N — > oo, £ finite, depends on N — |, the same combination 
as in the case of the harmonic oscillator, as shown by one of us and C. Quigg [fn]| . 

An exhaustive study of the various cases has been made by Feldman, Fulton, 
and Devoto [fHfl . We shall not give in the text the refined formulae obtained by 
these authors. (A further refinement is given in Appendix B.) We shall content 
ourselves with the leading terms. Here we express energies in terms of the number 
of nodes n = n r — 1 = N — £ — 1 and the angular momentum £. We have, defining 
E(N,£) = E(n,£), the following cases: 

i) £ large, n finite, —2<u<oo: 



E(n, £) ~ e(u) ^ 



1+ 2 



x 



1 



ii a) n large, £ finite pjj, < v < oo: 



+ g + ( n + 2 ' 



2v 
v+2 



(12) 



E(n, 



2^ l2 + ^ 



r(i + i)l n+ 2 I 



2v 
v+2 



(13) 



ii b) n large, £ finite /fjj/, — 2 < u < 0: 
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Figure 1: Energy levels in the Schrodinger equation (|TT|) for several power poten- 
tials. Levels are labeled as described in the text. 
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These formulae exhibit very smooth behavior of the energies as functions of n 
and £ and hence N and £. Further support for the smoothness in £ for fixed n is 
given in Appendix B, because [E(n,£)]» is a "Herglotz" function of (£ + 1/2) 2 
for v > 0. 

In the case i), for £ — > oo, n finite, one can write a systematic expansion of the 
energies in inverse powers of £ + | [13|. The first terms are 



1v 
v+2 



X 1 + 



2z/ n + 



i 



+ 



v{y-2) 



^17+2 £ + \ 12(z/ + 2)(£ + ±) 2 

« r i/ - 2 

+e> 



i 

+ 2 



ll-,)(n + -) 2 + ^- 



(£+±) 3 



(15) 



Notice that this expression becomes exact for v = 2. In Appendix B we show 
that (plj) can be used to improve the Feldman-Fulton-Devoto []12 asymptotic 
expression for large £ in a rather impressive way. Eq. (|15|) yields the behavior of 
Rn/ for £ — * oo, N — £ finite: 



R 



Nl 



{y-2f 4-3 v / ^T2-z/ 
6(2£ + l) ( Z y + 2)( v / ^7T2 + 2") 

2) 2 



(2£+l) 2 ■ (16) 

This expression shows that is indeed very close to unity for — 1 < v < 2, as 
long as £ > 2. 

Other limiting cases which can be treated analytically are v — > — 1 and v — > 2, 
for arbitrary n and £. 

In the case v — > — 1 it is sufficient to look at a perturbation of the Coulomb 
potential by a potential log(r)/r, since 



lim 



log(r) 



-i v + 1 r 
Calculations presented in Appendix C give 

N + £+ 1 



(17) 



i? 



iV + 



for i/ — > — 1. 

In the case v — > 2, one should, similarly, take a perturbation V = r 2 logr. 
Then one gets the surprise that the terms of order v — 2 vanish and 



i? 



1 + -2) 



(19) 
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At this point, we see that it may be advantageous to replace Eq. ( fTBD by 



(u-2) 2 4-3 v / ^T2-z/ 
RN ^ 1 + 6(N + l) (, + 2)0 + 2 + 2) ' (20) 



which preserves the large-£ asymptotic behavior and agrees with (jig) for v — > — 1 
and (0) for u -> 2. 

All this fits perfectly into the numerical calculations presented in Figs. 2 and 
3 for (N = 3, £ = 1), (iV = 4, £ = 1), and (JV = 4, £ = 2). It is striking that the 
asymptotic formula (|I~6| ) agrees quite well with the case in which n is smallest and 
£ largest, i.e., N = 4, £ = 2. In particular the asymptotic formula reproduces the 
fact that i? — 1 is negative for ~ 0.1 < v < 2. 

For completeness let us indicate that for v — > —2, R^j tends to +oo. Notice 
first that, by scaling, R^,e is independent of the strength of the power potential. 
We can always take V = —gr v with (£+ \) 2 > g > (£ — \) 2 . Then E N> i and E N j + x 
approach zero as v — ■> —2, while E^/-\ goes to — oo. This is because, for a state 
of angular momentum £, the limit Hamiltonian is, for v — > — 2, 

Ham d2 1 | _ (21) 

"" dr 2 4r 2 r 2 

The operator —(d 2 /dr 2 ) — l/(4r 2 ) is known to be positive, i.e., its expectation 
value for any reduced wave function vanishing at the origin is positive. Hence, if 
g < (£ + \) 2 the eigenvalues of Hi im are positive. However, for > v > —2, the 
Hamiltonian has negative eigenvalues (in fact infinitely many!) and therefore the 
lowest eigenvalue has to go to zero for v — > —2. On the other hand for g > (£+ 1) 2 
it is known that H\{ m is not lower-bounded, as can be seen by taking its expectation 
value with, for instance, a trial function u = for r < R m , (r — R m )2 +e (R M — r) for 
R m < r < R M , for r > R M , and letting R m and e go to zero. This not only proves 
that the ground state energy goes to — oo, but all radial excitation energies do so as 
well. If E n is a finite limiting value of the energy level of the nth radial excitation, 
there is a sequence of energies and wave functions (defined by integration of the 
Schrodinger equation from infinity) which approaches the limit energy and wave 
function. Let be the first nonzero limit of a node for v — > —2. Then taking 
r fe _x < R m < R M < r k , the energy corresponding to the interval (R m , Rm) with 
Dirichlet boundary conditions is higher than the one corresponding to (?"£._!, r^.) 
and goes to — oo, which is what we wanted to prove. If all nodes approach zero 
(which is in fact the case!) the proof obviously works, taking R m > r n . 

The opposite extreme case is v — > oo, i.e., the square well. Then the solutions 
of the Schrodinger equation are Bessel functions which can be approximated near 
their turning point by Airy functions. One gets, for instance, 



lim R e+2 ,e = fn f^T = 0.818709 ... , (22) 
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Figure 2: Ratios of spacings of levels with the same principal quantum number 
in power-law potentials V(r) ~ r v as a function of v. Solid curve: [E(3S) — 
E(2P)]/[E(2P) - E{1D)\\ Dashed: [E{AS) - E(3P)]/[E(3P) - E(2D)\; Dotted: 
[E(3P) - E(2D)]/[E(2D) - E(1F)\. Here levels are labeled by n r (S, P, D, F, . . .), 
where n r = N — £ is the radial quantum number. 
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Figure 3: Ratios of spacings of levels with the same principal quantum number in 
power-law potentials V(r) ~ r u near v = 2. Solid line: [E(3S) - E(2P)]/[E(2P) - 
E(ID)\] dashed line: [E(4S) - E(3P)]/[E(3P) - E(2D)]; dotted line: [E(3P) - 
E(2D)]/[E(2D) - E(1F)\. Levels are labeled as in Fig. 2. 
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where the at are successive zeros of the Airy function. 

IV. APPLICATIONS TO QUARKONIUM LEVELS 

Realistic quarkonium potentials appear to have a power-law behavior V ~ r u , 
@], with v not far from zero. For -0.1 < v < 0.1, R 3jl is 1.007 ± 0.002. Con- 
sequently, one can anticipate that energy levels with the same principal quantum 
number N will depend nearly linearly on the orbital angular momentum i. We ap- 
ply this result to two examples, one in charmonium and one in the upsilon family. 



A. Charmonium 2P states 

The IP levels of charmonium were identified nearly 20 years ago in electric 
dipole transitions from the ip'(2S) [|i"4jl . Recently they have been studied with 
great precision in proton-antiproton formation experiments Since they lie 

below threshold for decay to a pair of charmed mesons, they are quite narrow, 
facilitating their observation. The 3 Pi and 3 P2 levels, in particular, have substantial 
branching ratios to 7 + J/ip. 

Recent interest has focused on the possibility that one or more charmonium 2P 
states may be narrow enough to have a substantial branching ratio to 7 + ip'(2S) 
0, [5|]. This suggestion is motivated by a production rate for ip'(2S) in high- 
energy proton-antiproton collisions |16| which appears too high to be explained by 



conventional QCD mechanisms. It then becomes of some interest to predict exactly 
where the 2P levels should lie. If the 3 P2 level lies sufficiently low, its decay to 
DD, though kinematically allowed, will be suppressed by a large 1 = 2 centrifugal 
barrier, so that its branching ratio to 7 + ip'(2S) could be non-negligible. The 3 Pi 
level cannot decay to DD; it will be narrow if it lies below the DD threshold 
at 3.87 GeV. Reliable anticipation of the positions of the 2P levels may permit 
their discovery and study in the same low-energy direct-channel experiments which 
studied the IP levels so successfully [|15| . 

In order to anticipate the spin-weighted average 2P mass in charmonium, we 
need similar quantities for the 3S and ID levels. The masses of the 3 3 Si and l 3 Di 
states are quoted as 4040 ± 10 and 3769.9 ± 2.5 MeV, respectively Assuming 
hyperfine splittings in the S'-wave levels to scale roughly as 1/N [[|, and using the 
observed splitting of about 118 MeV between the 1 3 S\ and l 1 ^ charmonium levels, 
we expect the spin- weighted average 3S mass to be about 4030 ±10 MeV. We must 
rely on a specific calculation [^] of fine structure to estimate the spin-weighted 



average ID mass; the result is 3820 MeV. Thus we predict the spin-weighted 2P 
mass to lie around 3925 MeV. The prediction of a 2P level near the average of the 
3S and ID levels should be a feature of any smooth potential which reproduces 
charmonium and upsilon levels. 
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B. ID states of the upsilon family 



Electromagnetic transitions from the T(3S) to the x' b (2P) levels, followed by 
transitions to the Y(LD) levels and their subsequent decays, can give rise to char- 
acteristic photon spectra to which the CLEO-II detector is uniquely sensitive. 
It is also possible for the CLEO-II detector to scan in energy for the T(l 3 Di) level, 
whose leptonic width is expected to be one or two electron volts |LS]. We expect 
the spin- weighted T(1D) mass to be the same distance below the spin- weighted 2P 
mass (10.260±0.001 GeV) as the distance between this mass and the spin-weighted 
3S mass. Here we do not know the hyperfme splitting between the observed T(3S) 
at 10.3553 ± 0.0005 GeV and its spin-singlet partner. On the basis of the range 
of predictions for the IS* level and our assumption that hyperfme splittings scale 
as 1/N, we estimate the hyperfme splitting in the 35' system to be 20 MeV, give 
or take a factor of 2, and hence the spin-averaged 3S level to lie at 10.349 ± 0.007 
GeV. Thus, we would expect the spin-averaged T(1.D) level to lie around 10.17 
GeV. 

C. 2D and 3P states of the upsilon family 

It may be possible to detect 2D upsilon levels through direct-channel e + e~ 
annihilation or hadronic production ||. Potential models predict the spin-average 
of the 2D levels to lie around 10.44 GeV. Let us imagine that such a level has 
indeed been seen. Then, since the T(4S I ) state has a mass of about 10.58 GeV, 
the 3P level should lie midway between the two in a power-law potential, at 10.51 
GeV. Since this is below BB threshold, the 3P level should be narrow, as has been 
noted previously ||. This level, more easily produced in hadronic interactions than 
the T(3S), should then be able to populate the 3S level through electric dipole 
transitions, which may be detectable. Whichever level, 2D or 3P, is seen first, our 
equal-spacing rule permits the other's mass to be anticipated. 

D. Possible sources of error in predictions 

Explicit potential models predict a slightly lower T(1D) mass (ranging from 
10.15 to 10.16 GeV) than our prediction of 10.17 GeV, as a consequence of the 
change in the effective power of the potential with distance. The 3S — 2P splitting 
is more sensitive to the short-distance (more Coulomb-like) part of the potential, 
while the 2P— ID splitting is sensitive to longer-range effects, where the potential 
is expected to be closer to linear. For example, a typical Coulomb-plus-linear 
potential, such as V(r) = — 0.4/r + 0.16r, where r is in GeV -1 , and mj = 5 GeV, 
implies E(3S) — E(2P) around 15% less than E(2P) — E(1D). To a lesser degree, 
such changes in effective power of the potential may be visible in deviations of the 
charmonium 2P and upsilon 3P masses from our predictions. 

Coupled-channel effects can significantly distort predictions of potential models 
near pair production thresholds [ljj. We expect such effects to be most significant 
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in the anticipation of the charmonium 2P and upsilon 3P levels mentioned above, 
and less important for the upsilon ID level. 



V. SUMMARY 

Interesting questions remain to be settled in quarkonium systems. The nonrel- 
ativistic Schrodinger equation continues to provide a useful guide to the properties 
of such systems, with power- law potentials V ~ r u , v 0, permitting rapid an- 
ticipation of charmonium and upsilon properties. We have shown that in a wide 
range of power-law potentials, energy levels characterized by the same principal 
quantum number are approximately linear in orbital angular momentum £, with a 
coefficient which is negative as long as v > — 1. We have exhibited this behavior 
numerically, discussed the limiting cases of perturbations around the Coulomb and 
harmonic oscillator potentials, presented semiclassical results for the coefficient of 
£, and applied the results to the anticipation of several charmonium and upsilon 
levels which have yet to be observed. 
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APPENDIX A: PROOF OF THE BAUMGARTNER THEOREM 
AND EFFECT OF LINEAR PERTURBATIONS 

The Schrodinger equation for a pure Coulomb potential is taken to be 

To prove the Baumgartner theorem, we let 

&N,e — E<N,e — E Njt+1 . (Al) 

Then we have 

A N/ = J v (u 2 Ne - u 2 N/+1 ) dr , (A2) 

where the w's are Coulomb wave functions, un/ and un/+i are linked by raising 
and lowering operators: 

A± , d £ + 1 




±- — + 



dr r 2(£+l) 
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A e u N: i+i = aeu N>i , Afu Nj £ = a e u N: e +1 , (A3) 
with Aaj = (£ + 1)~ 2 — A^~ 2 . The the Coulomb Hamiltonian can be written as 

H e = AjAf- A{g+1)2 ■ 
Using these operators, one can transform A^t into 

/•OO 

a e A Ni = / V'u N>e+ iu N/ dr {AA) 
Jo 



and 



PCX) rOO j 

a^+iAjv^+i = / V'u N/+1 u N/+2 dr = / V^'u^m-i Aj +1 u N/+1 dr 

Jo Jo Ctl+l 



Hence, using the expression for At +1 and integrating by parts, 



a e+1 A N/+ i 



= f 

Jo 



2 r 2{£ + 2) 



Similarly 



and 



roc 

a e A N>e = / V 
Jo 



c? c A 



- f 
N/ ~ Jo 



+ 



V" t + l irl 



-A f u N . 



un,£+i 



V' + 



2(e + i) 



2 r 

Combine now equations (A5) and (A6): 

{£ + l)a 2 e A N/ -(£ + 2)a 2 e+1 A N/+1 

2£ + 3 



u2 N,t+idr 



= J ((2€ + 3)^ + (2€ + 3)^) 



(A5) 



(A6) 



J AVu% /+1 dr (A7) 



which is what we wanted to prove. 

Now we can use Eq. (A6) to calculate Ajy/ for the case of a linear perturbation, 
Ar, to the Coulomb potential. We get 



ajA Nt i = A 
and, using the virial theorem, 



1 1 



r 2(£+l) 



a e A Nti = -A 



2(£ + l)\E^jlT h \ ~ 



2(£+l) 



A N 2 - (£+lf 
2 N 2 {£+1) 



(A8) 
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Hence, from the definition of a£, E^ — E^^+i = 2A(£ + 1), and, therefore, for 
V = -(l/r) + Xr, 

E3S — E 2 p \ 1 



lim ^ ^ = - . (A9) 



APPENDIX B: SOME COMMENTS ABOUT 
POWER POTENTIALS IN THE LARGE £ LIMIT 



i) Equation flTq ) has been obtained in Ref. fl3|| , but in this reference, the 

— 1 /2 

expansion parameter was [£{£ + 1)] • The advantage of re-expressing things in 
an expansion in (£+\)~ 1 is that for v = 2 the expansion stops and gives the exact 
answer. 



ii) Equation fll2 ) gives only the leading term of the large £ behavior obtained by 
the semi-classical approximation, while Feldman, Fulton and Devoto (FFD) give 



E(n, 



~ £ (V 




v+2 



1 



1 , 



- + (n+-)vV + 2 



2v 
v+2 



X 



24 



+ 1 + 



n 



v 



2v 
v+2 



(Bl) 



When one expands this expression in inverse powers of (£ + =) 1 and compares 



with Eq. (15) which is the beginning of a systematic asymptotic expansion in 



(£ + |) 1 one finds a small difference which can be corrected by subtracting 1/12 
from (n+ |) 2 in the second bracket: 



E(n, 




v+2 



v 

1+ 2 



1 

+ 2 + 



n + -)Vi> + 2 



2v 
v+2 



X 



1 - 



12 



24 



(52) 



It turns out that the small change made in the second bracket makes this new 
expression extremely accurate especially for small n and in particular for n — 0. 
Then for v = —0.5, v = 0.5, v — 1 we find by numerical tests that the relative 
error of Eq. (B2) is less than 10 -4 for n — 0, arbitrary £ > 0. For n — 1 it is 
less than 5 x 10 -3 . We believe that this holds on the whole interval — 1 < v < 2. 
Examples of the accuracy of Eq. (B2) for v = —0.5 and v = 1 are shown in Tables 
1 and 2. 
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Table 1: Comparison of exact and approximate [Eq. (B2)] energy levels for the 
Schrodinger equation (Q) in a potential V(r) = — r~ 1//2 . 



Exact 


N = 


- 1 


N = 2 


N = 3 


N = 4 


S 


-0.43? 


3041 


-0.263203 


-0.197558 


-0.161705 


P 






-0.286611 


-0.209800 


-0.169416 


D 








-0.221506 


-0.176817 


F 










-0.184005 


Eq. (B2) 


N = 


= 1 


N = 2 


N = 3 


N = 4 


S 


-0.43? 


^043 


-0.264647 


-0.199228 


-0.163352 


P 






-0.286615 


-0.210156 


-0.170025 


D 








-0.221507 


-0.176947 


F 










-0.184006 



Table 2: Comparison of exact and approximate [Eq. (B2)] energy levels for the 
Schrodinger equation (|) in a potential V(r) = r. 



Exact 


N -- 


= 1 


N = 2 


N = 3 


iV = 4 


S 


2.33? 


3107 


4.087949 


5.520560 


6.786708 


P 






3.361255 


4.884452 


6.207623 


D 








4.248182 


5.629708 


F 










5.050926 


Eq. (B2) 


N -- 


= 1 


N = 2 


N = 3 


N = 4 


S 


2.33? 


3231 


4.066542 


5.479175 


6.727770 


P 






3.361231 


4.874358 


6.183454 


D 








4.248160 


5.623973 


F 










5.050910 
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iii) The smoothness of E(n,£) for fixed n, as a function of £ can be connected 
to the fact that 



E(n,£) {V+2) ' U is a "Herglotz" function of {£ + \) 2 . A Herglotz 



function is defined by 

z f° lmH(z')dz' 

z'(z' - z) 



H{z) = A + Bz + - [ 

7T J-, 



with ImH(z') > 0. It has the characteristic property lmH(z)/lmz > 0, and cannot 
grow faster than z in any complex direction. Furthermore, it is easy to see that it 
is concave, i.e. d 2 H(z)/dz 2 < 0. 

It has been shown by one of us (A. M.) and Harald Grosse [20] that E(n,l) for 
V = r v , v > is analytic in the variable z = {£ + |) 2 in a cut plane, with the cut 
running from z = — oo to z — 0. For positive z, E is real. Furthermore for Imz > 
we have < Arg E < ^vr, hence < Arg (£)(^+ 2 )/^ < n . Hence E (,/+2 ^ u is a 
Herglotz function of z, and, as we said, is, in particular, concave. Let us illustrate 
the usefulness of this remark by considering V = r 4 . We have, for instance 

£(4,0) = 44.0 £(4,1) = 50.1 £(4, 2) =56.4 . 
Using the concavity of E z l 2 in {£ + |) 2 we get from E(4, 0) and E(A, 2) the bound 

£(4,1)>48.3 . 



APPENDIX C: 
BEHAVIOR OF R N/ NEAR v = -1 AND v = 2 

To study the behavior of Rn,£ in the neighborhood of v = — 1 we have to take 
a perturbation of the Coulomb potential which is v — (log r)/r, but we shall first 
look at a perturbing potential which is v — logr. Then according to equation (A6) 
of Appendix A we have 

2 a / 1 „ £ + 1 , 1 A 2£ + 3 / 1\ 1 /1\ 

^ = \2 V ~ — V + 2^^ ) N/+1 = ~— \- 2 ) N ^ + WTT) \r) N/+1 

= l¥ ln fixcd + 2(£ + 1) ^ 2E) = APP {l+l ~ n) ■ 
Therefore, for v = logr, 

A « = N WTTTT) = L t«h-«h»)logrdr . (C2) 
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Now, following a standard strategy of multiplying the Schrodinger equation by 
■ulogr and u'rlogr and integrating, it is possible to obtain the following identity: 



Hence combining (C2) and (C3) we get 

logr 



i u N,e u N,e+i 



dr = 



and hence 



R 



2N 2 (N + £+l) 
Em^ - E N , N + £ + l 



N,e 



En,£ — Enp + i 



N + 



(C4) 



(C5) 



Now we turn to the neighborhood of v — 2. The harmonic oscillator Hamilto- 
nian is taken to be 



d 2 £{£+!) r 2 



(C6) 



dr 2 ' r 2 4 

Then, fori/ = r 2 /4, we have E(n,£) = 2n+i+\, and E(n + 1,£-1) -E(n,£) = I. 
We consider now the perturbation 



With 



we have 



St = J q v (u 2 n/ - dr 



R ±- + d t + lr 



Bg~U n+ ij — 7 n Wn,m , B t — 7„, U n+ ij , 



(C7) 



(C8) 



(C9) 



with 7„ = 2n + 2. Then, by repeated use of these raising and lowering operators 
one can get 



*-*-n[( 



2 \ ' 2 2 

rM \ U In t\ 2 r 2 



7n^-l = / - 

Jo r 

where u = u n ^. Then we get 



2 \ ' 2 2 

ru \ u „ 9 r 9 
H £w H M 



dr , 



5/ — 5, 



2n + 1 2 



4n(n + 1) 



ru dr 



2n(n + 1) Jo r 



-E(n,. 



u dr 



(CIO) 
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For v = log r and for v = r 2 log r all these integrals can be calculated and one 
finds in both cases 5e — <^_i = 0. Hence if the total potential is 

V = r v = r 2 + r v - r 2 ~ r 2 + {v - 2)r 2 logr , 

then 5e — Se-i is of higher order, i.e. Oiv — 2) 2 , which fits with the numerical 
observations of Fig. 3. 

References 

[1] D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 (1973); Phys. Rev. 
D 8, 3633 (1973); 9, 980 (1974); H. D. Politzer, Phys. Rev. Lett. 30, 1346 
(1974); Phys. Rep. 14C, 129 (1974). 

[2] Y. Nambu, Phys. Rev. D 10, 4262 (1974) 

[3] M. Machacek and Y. Tomozawa, Ann. Phys. (NY.) 110, 407 (1978); C. Quigg 
and J. L. Rosner, Phys. Lett. 71B, 153 (1977); C. Quigg and J. L. Rosner, 
Comments on Nucl. Part. Phys. 8, 11 (1978); C. Quigg and J. Rosner, Phys. 
Rep. 56C, 167 (1979); A. Martin, Phys. Lett. 93B, 338 (1980); 100B, 511 
(1981); C. Quigg, in Proceedings of the 1979 International Symposium on 
Lepton and Photon Interactions at High Energies, Fermilab, August 23-29, 
1979, ed. by T. B. W. Kirk and H. D. I. Abarbanel (Fermi National Accel- 
erator Laboratory, Batavia, IL, 1979, p. 239; A. Martin, in Heavy Flavors 
and High Energy Collisions in the 1-100 TeV Range, edited by A. Ali and L. 
Cifarelli (Plenum, New York, 1989), p. 141; A. Martin, CERN report CERN- 
TH. 6933/93, 1993, invited lecture at the 30th Course of the International 
School of Subnuclear Physics, Erice, July, 1992; A. K. Grant, J. L. Rosner, 
and E. Rynes, Phys. Rev. D 47, 1981 (1993). 

[4] P. Cho, M. B. Wise, and S. P. Trivedi, Phys. Rev. D 51, 2039 (1995). 

[5] F. E. Close, Phys. Lett. B 342, 369 (1995). 

[6] W. Kwong and J. L. Rosner, Phys. Rev. D 38, 279 (1988). 

[7] See, e.g., W. Buchmuller and S.-H. H. Tye, Phys. Rev. D 23, (1981); R. 
McClary and N. Byers, Phys. Rev. D 28, 1692 (1983); S. N. Gupta, S. F. 
Radford, and W. W. Repko, Phys. Rev. D 26, 3305 (1982); 28, 1716 (1983); 
31, 160 (1985); 34, 201 (1986); S. N. Gupta, C. J. Suchyta, and W. W. Repko, 
Phys. Rev. D 39, 974 (1989); W. Buchmuller and S. Cooper, in High Energy 
Electron- Positron Physics, edited by A. Ali and P. Soding (World Scientific, 
Singapore, 1988), p. 488; L. P. Fulcher, Phys. Rev. D 37, 1258 (1988); 42, 
2337 (1990); 44, 2079 (1991); L. P. Fulcher, Z. Chen, and C. K. Yeong, Phys. 
Rev. D 47, 4122 (1993); D. Lichtenberg et al, Zeit. Phys. C 41, 615 (1989). 



19 



[8] H. Grosse and A. Martin, Phys. Lett. 134B, 368 (1984); B. Baumgartner, H. 
Grosse, and A. Martin, Phys. Lett. 146B, 363 (1984); Nucl. Phys. B254, 528 
(1985). 

[9] B. Baumgartner, Ann. Phys. (N.Y.) 168, 484 (1986). 

[10] A. Grant and J. L. Rosner, Am. J. Phys. 62, 310 (1994). 

[11] C. Quigg and J. L. Rosner, Phys. Rep. 56, 167 (1979). 

[12] G. Feldman, T. Fulton, and A. Devoto, Nucl. Phys. B154, 441 (1979). 

[13] H. Grosse and A. Martin, Nucl. Phys. B254, 528 (1985). A misprint in the ex- 
pression for energy levels (a factor 1/24) is corrected by A. Martin, in Recent 
Developments in Mathematical Physics (Proceedings of the XXVI Interna- 
tionale Universitatswochen fur Kernphysik, Schladming, Austria, Feb. 17-27, 
1987), edited by H. Mitter and L. Pittner (Springer Verlag, New York, 1987), 
p. 53. 

[14] W. Braunschweig et al, Phys. Lett. 57B, 407 (1975); G. J. Feldman et al, 
Phys. Rev. Lett. 35, 821 (1975); W. Tanenbaum et al, Phys. Rev. Lett. 35, 
1323 (1975), Phys. Rev. D 17, 1731 (1978) and references therein. 

[15] Fermilab E760 Collaboration, A. Ceccuci et al, Nucl. Phys. A558, 259c 
(1993); T. A. Armstrong et al, Phys. Rev. D 48, 3037 (1993); Phys. Rev. 
Lett. 70, 2983 (1993). 

[16] CDF Collaboration, T. Daniels et al, Fermilab report FERMILAB-CONF- 
94/136-E, 1994 (unpublished); CDF Collaboration, K. Byrum et al, Fermilab 
report FERMILAB-CONF-94/325-E, Sept., 1994, presented at 27th Interna- 
tional Conference on High Energy Physics, Glasgow, Scotland, July 20 - 27, 
1994. 

[17] Particle Data Group, L. Montanet et al, Phys. Rev. D 50, 1173 (1994). 

[18] P. J. Moxhay and J. L. Rosner, Phys. Rev. D 28, 1132 (1983); P. J. Moxhay, 
Ph. D. Thesis, University of Minnesota, 1983 (unpublished). 

[19] N. Byers and E. Eichten, Phys. Rev. D 42, 3885 (1990); in EPS - High 
Energy Physics '89 (Proceedings of the International Europhysics Conference 
on High Energy and Particle Physics, Madrid, Spain, Sept. 6 - 13, 1989), 
edited by F. Barreiro et al (North-Holland, Amsterdam, 1990), Nucl. Phys. B 
Proc. Suppl. 16, 281 (1990). 

[20] H. Grosse and A. Martin, Phys. Rep. 60, 341 (1980). 



20 



